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INTRODUCTION 


Gravimeters  are  currently  being  used  on  shipboard  surveys  to  measure  the 
earth's  gravitational  field.  The  gravimeter  measures  vertical  accelerations  due  to 
nongravitational  forces  acting  on  the  vehicle,  since  mean  sea  level  offers  an 
excellent  vertical  reference  and  the  ship  is  moving  at  a  low  velocity,  these  ac¬ 
celerations  are  averaged  over  a  time  interval  to  produce  a  value  very  nearly  equal 
to  the  negative  of  the  gravitational  acceleration  magnitude.  For  ships  moving  with 
a  speed  of  about  10  kn,  the  error  due  to  this  averaging  out  of  vertical  motion  is 
estimated  to  be  a  few  tenths  of  a  mgal  or  less  (Reference  1).  This  is  small  when 
compared  to  overall  shipboard  survey  accuracy  of  about  2  mgal  (Reference  1). 

For  a  higher  velocity  survey  vehicle  (e.g.,  an  airplane),  the  averaging  time 
becomes  much  smaller  and  is  the  limit  factor  for  survey  accuracy  (Reference  2).  It 
becomes  necessary  to  use  an  external  vertical  reference.  Here,  the  gravimeter  data 
are  subtracted  from  the  vertical  reference  and  then  estimated.  Up  to  now,  suf¬ 
ficiently  accurate  vertical  referencing  has  been  difficult  to  obtain,  especially 
over  land  (Reference  2). 

This  report  investigates  the  use  of  a  Global  Positioning  System  (GPS)  receiver 
for  vertical  referencing.  (See  Reference  3  for  GPS  background.)  This  investigation 
is  motivated  by  recent  advances  in  GPS  receiver  capability  and,  specifically,  the 
development  of  the  GEOSTAR*  receiver  (References  4  and  5).  This  receiver  tracks  up 
to  four  satellites  simultaneously  and  offers  the  potential  for  very  accurate  verti¬ 
cal  measurements.  Simulated  and  actual  measurement  accuracies  of  the  GEOSTAR  are 
given  in  the  next  section. 

The  purpose  of  this  study  is  to  determine  if  the  currently  available  GPS  re¬ 
ceivers,  specifically  the  GEOSTAR  receivers,  are  sufficiently  accurate  to  be  used 
for  vertical  referencing  in  gravity  surveys.  Only  the  error  in  vertical  referencing 
due  to  the  receiver  itself  is  used  in  the  study.  Consequently,  it  is  assumed  that 
all  values  used  are  known  perfectly  except  for  the  receiver  phase  measurement  noise. 
This  assumption  enables  values  of  nongravitational  vertical  position,  velocity,  or 
acceleration  determined  from  the  gravimeter  to  be  subtracted  directly  from  corres¬ 
ponding  GPS  determined  values.  The  residuals  of  these  subtractions  are  due  only  to 
the  gravity  effect  and  measurement  noise  of  the  receiver.  Therefore,  the  particular 
vehicle  motion  does  not  need  to  be  modeled  statistically,  instead,  only  the  gravity 
anomalies  and  receiver  measurement  noise  need  to  be  statistically  modeled.  The 
remaining  accuracy  analysis  assumptions,  the  statistical  model  for  the  gravity 

•Texas  instruments,  Inc.,  is  designing  and  building  the  GEOSTAR  receiver  under  con¬ 
tract  to  the  Applied  Research  Laboratories  (ARL).  NSWC  is  developing  the  receiver's 
control  and  navigation  software.  DMA  National  Oceanic  &  Atmospheric  Administration 
(NOAA)  and  U.S.  Geological  Survey  (USGS)  are  funding  the  development. 
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anomalies,  and  the  relationships  for  the  optimum  smoothing  filter  and  the  corres¬ 
ponding  mean  square  estimation  error  are  discussed  herein. 


RECEIVER  MEASUREMENT  ACCURACY 


The  GEOSTAR  receiver  was  designed  to  make  very  precise  phase  measurements  of 
the  GPS  satellite  transmitted  signals.  Simulations  have  been  performed  at  Texas 
Instruments  (Reference  6)  to  predict  the  expected  phase  measurement  accuracy  for 
several  dynamic  environments.  These  values  are  listed  in  Table  1  for  stationary, 
low  dynamic  (boat  or  jeep),  and  medium  dynamic  (helicopter  or  airplane)  cases.  For 
all  of  these  cases,  the  maximum  expected  receiver  phase  measurement  error  standard 
deviation  is  about  1  cm.  These  values  were  obtained  using  a  simulation  model  and  do 
not  include  all  real  world  effects  (e.g.,  vehicle  vibrations).  A  phase  measurement 
is  equivalent  to  the  range  to  the  satellite  with  the  addition  of  an  unknown  bias. 
However,  since  derivatives  or  changes  in  range  are  of  interest,  these  bias  values 
are  unimportant. 


TABLE  1 .  EXPECTED  GEOSTAR  MEASUREMENT  ACCURACY* 

Phase  Vertical** 

Maximum  Measurement  Position 
Velocity  Error  Error 

_ Case _  (m/s )  (m)  (m) 

Stationary  0  0.0015  0.003 

Low  Dynamic  25  0.003-0.005  0.006-010 

(Boat  or  Jeep) 

Medium  Dynamic  200  0.003-0.01  0.016-0.02 

(Helicopter  or 

Airplane) 

♦Values  quoted  for  accuracy  are  theoretical  instru¬ 
mentation  values.  Subsequent  hardware  tests  have 
verified  that  these  values  are  real  within  a  factor 
of  two  or  three. 

**RMS  Vertical  Dilution  of  Precision  (VDOP)  =  2 


More  specifically,  the  vertical  accelerations  are  to  be  estimated  from  the 
receiver  biased  range  measurements.  The  accuracy  with  which  the  vertical  component 
can  be  computed  depends  on  the  number  of  satellites  in  view  and  their  relative  loca¬ 
tion  with  respect  to  the  receiver  antenna.  The  decrease  in  accuracy  due  to  the 
locations  of  the  satellites  in  view  is  called  Vertical  Dilution  of  Precision  (VDOP). 
A  value  of  VDOP  *  2,  chosen  from  Reference  3,  is  the  50th  percentile  value  and  close 
to  the  root  mean  square  value.  This  value  is  based  on  a  range  error  budget  of  3.6 
to  6.3  m,  a  24-satellite  baseline  constellation,  and  a  five-degree  satellite  mask 
angle.  Note  that  although  the  assumed  constellation  of  Reference  3  is  not  the 
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currently  planned  1 8-satellite/6-plane  constellation,  the  chosen  value  of  VDOP  is 
still  thought  to  be  reasonable.  The  GEOSTAR  receiver  also  measures  pseudorange.  It 
is  expected  that  the  resulting  range  measurements,  after  atmospheric  corrections, 
would  be  well  within  the  above  assumed  range  error  budget.  Multiplication  of  the 
VDOP  value  by  the  range  instrumentation  accuracy  produced  the  estimated  vertical 
measurement  accuracy.  The  vertical  measurement  accuracies  are  also  given  in  Table  1 

Measured  average  velocity  values  of  the  prototype  GEOSTAR  are  presented  in 
Reference  4.  These  values  are  taken  from  a  stationary  antenna,  but  the  tracking 
filter  parameters  were  not  set  for  the  stationary  case.  The  phase  measurements  were 
differenced  at  4-s  intervals  to  form  the  average  velocity  values.  This  was  done  for 
the  four  GPS  satellites  available  on  both  the  L  and  L  channels.  The  average 
standard  deviation  of  the  relative  velocity  measurement  error  computed  from  these 
data  was  0.003  m/s.  This  value  is  over  a  relatively  short  data  interval  of  9  min, 
but  it  also  includes  all  other  measurement  errors  in  addition  to  the  receiver  mea¬ 
surement  errors.  A  corresponding  vertical  velocity  error  standard  deviation  is 
0.006  m/s. 

The  above  values  of  vertical  position  (with  a  bias)  and  vertical  velocity  mea¬ 
surement  accuracies  indicate  the  current  capabilities  of  a  GPS  receiver.  These 
values  will  be  used  later  in  the  accuracy  analysis  to  determine  gravity  anomaly 
measurement  accuracy. 


ACCURACY  ANALYSIS  ASSUMPTIONS  AND  DEFINITIONS 


The  following  assumptions  are  made  for  the  accuracy  analysis: 

1.  The  survey  vehicle  is  on,  or  near,  the  earth's  surface  and  is  traveling  in 
a  constant  direction  at  a  constant  speed.  Uiese  are  reasonable  assumptions, 
since  survey  vehicles  would  generally  travel  this  way. 

2.  This  analysis  is  to  determine  if  the  receiver  measurements  are  sufficiently 
accurate  for  a  gravity  anomaly  survey.  All  other  values  are  assumed  to  be 
known  perfectly.  The  receiver  vertical  measurement  error  is  statistically 
modeled  as  white  noise  (i.e.,  constant  power  over  all  frequencies). 

3.  The  GPS  vertical  direction  is  normal  to  the  earth's  ellipsoid.  The  gravi¬ 
meter  vertical  direction  is  normal  to  the  earth's  geoid.  This  small  angular 
difference,  normally  about  a  few  arc  seconds,  is  assumed  to  have  a  negligible 
effect  on  the  magnitudes  of  the  acceleration  values.  In  other  words,  the  mag¬ 
nitude  of  the  acceleration  in  one  vertical  direction  is  very  nearly  equal  to 
the  magnitude  of  the  component  of  acceleration  in  the  other  vertical  direction. 
Also,  perfect  modeling  of  the  vertical  coupling  between  the  gravimeter  and  the 
GPS  antenna  is  assumed. 

4.  The  data  are  assumed  to  be  sampled  sufficiently  fast  so  that  continuous 
signal  analysis  procedures  can  be  used.  Hie  sampling  interval  chosen  for  the 
survey  gravity  data,  discussed  in  the  next  section,  is  1  nmi  (Reference  7). 

The  Nyquest  frequency  (the  reciprocal  of  twice  this  sampling  interval) 
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represents  the  highest  frequency  of  interest  in  the  gravity  data  (Reference  8). 
A  sampling  interval  of  one  sample/s  in  a  vehicle  at  the  maximum  velocity  of  the 
study,  100  m/s,  would  have  a  cutoff  frequency  of  18  times  the  practical  gravity 
cutoff  frequency.  The  GEOSTAR  receiver  is  capable  of  taking  measurements  of  up 
to  25  times/s.  Therefore,  the  sampling  rate  appears  to  be  sufficient  for  this 
assumption. 

5.  The  gravity  anomalies  can  be  statistically  modeled.  Here,  stationarity  and 
other  assumptions  (Reference  9)  are  required.  Ihe  statistical  model  chosen  is 
described  later  in  this  report. 


The  following  definitions  are  used  in  the  analysis : 
g(t)  =  true  gravity  magnitude 

=  (total-vertical  acceleration  determined  from  GPS)-(nongravitational 
vertical  accelerations  measured  by  the  gravimeter) 


Ag(t)  =  gravity  anomaly 

=  g(t)-( ellipsoidal  model  gravity  magnitude) 

f(t)  =  vertical  position  due  to  g(t)  (with  bias) 

=  JT  Ag(X)dXdt 

z(t)  =  receiver  vertical  positioning  measurement  (white)  noise 


GRAVITY  ANOMALY  STATISTICAL  MODEL 


The  statistical  model  chosen  for  the  gravity  anomalies  is  given  by  the  auto- 
orrelation  function  (Reference  9) 

*AgAg(T)  =  °Ag  (1  +  0  M  “  M 2 }  exp  (  '  0  M  >  (1) 

here 


B  =  v/D 


v  =  constant  speed  of  the  survey  vehicle 
D  =  correlation  length 


°Ag 


=  gravity  anomaly  standard  deviation 


ote  that  since  the  vehicle  is  moving  in  a  constant  direction  and  at  a  constant 
peed,  the  radial  change  in  the  distance  parameter  of  the  autocorrelation  function 
f  Reference  9  is  converted  to  a  change  in  the  time  parameter.  This  analysis  is 
one  for  several  values  of  vehicle  speeds.  A  value  of  oAg  =  35.0  mgal  is  chosen  as 
worldwide  value  for  the  gravity  anomaly  standard  deviation  (Reference  10).  A 
alue  of  the  correlation  length  D  =  2.76x1  cf*  m  is  chosen,  since  it  matches  the  higher 
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OPTIMUM  SMOOTHING  FILTER 


The  problem  can  be  stated  as  follows:  Given  the  measured  vertical  position  due 

Agtt), 

z(t)  =  f(t)  +  n(t)  (3) 

:h  f(t)  and  n(t)  as  defined  above,  optimally  estimate 

Ag(t)  =  f"(t)  (4) 

:h  that  the  mean  square  error  of  the  estimate  is  minimized. 

The  measurement  noise  is  assumed  to  be  white.  Therefore, 

4>nu<0))  =  n2  (5) 

sre  n  is  a  constant  parameter  of  the  analysis.  From  Equation  (4), 

*ff(a))  =  ^AgAg^/w4  «»> 

ice  the  measurement  is  completely  uncorrelated  to  the  gravity  anomalies, 

4>fn(io)  =  *ne(u>)  =  0  (7) 

The  optimum  smoothing  filter  is  given  by  (see  Reference  12  and  13  for  details) 


H(ju>) 


=  »Aqz(oi) 

*zz<w) 


(8) 


re 


%z(a))  =  w2  $fz(w) 

*f,<«)  =  4>ff(w)  +  $fr)(w) 

*ZZ{10)  =  *ff(u)  +  *nn(w) 
sequertly, 


2a 


*A  (w)  = 

AcJz  m2 


2  3 

A„e 


r  p2  +  .^.2 1 

L  ( e2  +  uj2  > 3  J 

(9) 

v.‘.  ' 

zz 


2  a.  63 

(w)  = - Ja - 

w4 


6 2  ±  5U)2 
(82  +  w2)3 


+  n‘ 


(10) 
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:itution  of  Equations  (9)  and  (10)  into  (8)  gives 


(82  +  id2)3  82  +  5u>2  +  Ku>4  ( B2  +  u)2)3 


RESULTS 


The  parameter  values  previously  defined  are  substituted  into  Equation  (13)  and 
rically  integrated  for  different  values  of  noise  standard  deviations,  n,  and 
:le  speeds,  v.  The  results  are  shown  in  Figure  2.  The  estimation  error  is 
;ed  for  vehicle  speeds  of  1  (about  2  kn),  10,  20,  and  100  m/s.  As  the  vehicle's 
l  increases,  the  averaging  time  decreases  and  the  estimation  error  increases, 
jravity  anomaly  estimation  error  due  to  receiver  measurement  noise  is  less  than 
il  for  vertical  measurement  errors  of  5  cm  or  less  and  vehicle  speeds  up  to  100 
From  the  GEOSTAR  receiver  data  described  in  the  RECEIVER  MEASUREMENT  ACCURACY 
.on,  this  seems  to  be  within  the  receiver's  measurement  accuracy  capabilities. 


A  similar  plot  is  shown  in  Figure  3.  This 
rsxs  with  measurements  in  vertical  velocity 
;sponding  measurement  noise  is  now  in  units 
s  as  good.  However,  if  it  could  be  assumed 
nents  of  the  receiver  as  given  in  Reference 
.ty  anomalies  due  to  this  measurement  noise 
:le  velocities  of  100  m/s  or  less. 


is  the  result  of  repeating  the  above 
instead  of  vertical  position.  The 
of  m/s.  Here,  the  results  are  not 
that  the  noise  on  the  velocity  mea- 
4  is  white,  the  error  in  estimating 
would  be  less  than  2  mgal  for 


CONCLUSION 


The  results  shown  in  Figures  2  and  3  indicate  that  the  GEOSTAR  receiver  itself 
Ld  be  capable  of  making  sufficiently  accurate  measurements  to  be  used  as  a 
.cal  reference  on  gravity  anomaly  survey  vehicles  moving  at  speeds  of  100  m/s 
kn)  and  less.  Other  vertical  referencing  error  sources  will  be  discussed. 
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The  GPS  satellite  errors  contributing  to  survey  inaccuracy  include  orbit  error, 
antenna  vibration,  and  clock  error.  The  effect  of  these  errors  could  be  reduced  by 
having  additional  receivers  on  the  ground  at  precisely  known  positions  determined 
from  GPS  and  within  several  hundred  miles  of  the  survey  vehicle. 

Two  important  GPS  measurement  error  contributors  are  the  receiver  clock  error 
and  the  atmospheric  correction  errors.  On  the  short  term,  these  errors  seem  to  be 
sufficiently  small.  Over  time  spans  of  a  few  seconds,  the  change  in  the  clock  error 
and  the  change  in  the  atmospheric  error  are  very  small.  However,  the  longer  term 
changes  may  significantly  affect  survey  accuracy.  Also,  since  inertial  forces 
introduce  fluctuations  in  the  clock,  acceleration  compensation  in  the  oscillator  and 
special  shock  mounting  appear  to  be  necessary  (Reference  14). 

In  addition  to  the  above  error  sources,  accurate  compensation  for  movement  of 
the  electrical  phase  center  within  the  antenna  must  be  accomplished.  Also,  the 
coupling  between  the  GPS  antenna  and  the  gravimeter  must  be  modeled  very  well.  The 
errors  of  currently  available  gravimeters,  once  the  biases  are  removed,  appear  to  be 
small  compared  with  other  errors  in  the  survey  procedure  (Reference  1). 

In  addition  to  determining  vertical  acceleration,  the  GPS  receiver  is  capable 
of  determining  position.  Importantly,  horizontal  velocity  measurements  can  also  be 
obtained.  The  Horizontal  Dilution  of  Precision  (HDOP)  is  slightly  better  than  the 
VDOP  previously  discussed.  These  measurements  should  reduce  the  effect  of  the 
coriolis  error,  which  is  currently  the  limiting  factor  for  accuracy  of  the  shipboard 
surveys  (Reference  1). 
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